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We construct a dynamical model for the time evolution of the boson-fermion coexistence system. 
The dynamics of bosons and fermions are formulated with the time-dependent Gross-Pitaevsky 
equation and the Vlasov equation. We thus study the monopole oscillation in the bose-fermi mix- 
ture. We find that large damping exists for fermion oscillations in the mixed system even at zero 
temperature. 

PACS numbers: 32.80.Pj,67.57.Jj, 51.10. +y 

I. INTRODUCTION 

Over these last several years, there has been significant progress in the production of ultracold gases which 
t e aliz o Bo— co— <BEC S ) flflflQ, — e atomic Pe tmi gases fl, - B„ 8e - 
Fermi mixtures |6j. These systems offer great promise for studies of new, interesting driven-phenomcna. 
One of the most exciting themes in recent physics is to study the time-dependent dynamical motion 
of trapped atoms, such as collective oscillations [2( , quantum vortices |7( and atomic novae |8( ■ From a 
theoretical point of view, these phenomena are very important, as they allow us to construct and examine 
the transport theory in finite many body systems. Many body atomic systems are good probes for such 
study because the fundamental interactions are well understood and can be changed using the Feshbach 



resonance 



Of the above phenomena collective oscillations are very sensitive to properties of the system and therefore 
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are condensates for studying time-dependent dynamics. The collective oscillation of BECs have been 
studied experimentally 1£| and theoretically Furthermore, there have been theoretical studies of 

rmal phase ^| and in the superconducting phase ^]] , as well as the Bose- Fermi 



the Fermi gas in the normal 



mixed gases 



14, 
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17f and the BEC-BCS crossover 
Collective motions are usually studied with the random phase approximation (RPA). The RPA however, 
can only describe minimal vibrations around ground state; for example it is shown in Ref.|l5| that the 
radial variation is only about 0.05 %, while the amplitudes of actual experimental can be as large as 10 
%. 

To describe collective oscillations with larger amplitudes we need to calculate the time evolution of the 
system using the time-dependent Gross-Pitaevsky (TDGP) and time-dependent Hartree Fock (TDHF) 
equations. In atomic gas systems, however, the fermion number is greater than one thousand, but it is 
not easy to solve the TDHF equations with so many fermions and limited computer resources. For system 
with large fermion number, on the other hand, the semi-classical approximation (h — - > 0) becomes useful, 
and the Vlasov equation, corresponding to the semi-classical approximation of the TDHF equations, can 
be used. 

n 

In nuclear physics |19j the Vlasov-Uhling-Uhlcnbeck (VUU) approach, which is the Vlasov plus two 
body collision with the Pauli blocking, has been introduced, and very nicely explained many kinds of 



observables in heavy-ion collisions 



20]. In addition, this approach has been developed for the Lorentz 



covariant framework including the non-local mean-fields 



i ap 
2l|. 



As for the boson gas the time evolution of 



condensed bosons and thermal bosons have also been described with the time-dependent Gross-Pitaevsky 

n n 

equation and the Boltzmann equation |22j. Furthermore C. Menotti et al. [23j studied the expansion of 
a Fermi gas by calculating the Vlasov equation with the scaling method. 

In this work we construct a transport model including the condensed bosons and fermions based on the 
TDGP and Vlasov equations. As a first step we study the collective monopole motion of a the spherically 
symmetric system at zero temperature. This calculation is a good test for the comparison between the 
RPA and direct time evolution calculations. In addition one pays attention to the monopole oscillation 
from the point of the view of the collapse so as to know compressing processes of systems |2J| . 
In section 2, we describe our transport model. In Sec. 3 we give our results for monopole oscillations in 
the bose-fermi mixture, and discuss their properties. We then summarize in Sec. 4. 
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II. FORMALISM 

Here we briefly explain our formalism. In this work we consider the spherical boson-fermion mixing gas 
at zero temperature. 

First we define the Hamiltonian for boson-fermion coexistence systems as follows. 

H = H B +H F + H (1) 



with 



H B = I d 3 x 
Ilr - I ./ './• 



-\^{x)V 2 <p{x) + laV^M*) + ^{^(xM*)} 2 
1 -i/j^(x)V 2 ip(x) + ^-mfUjjx 2 il>\x)ip(x) 



(2) 
(3) 



Hbf = hsF J d 3 x [tf (x)<j)(x)^ {x)ip(x)] , (4) 
where <j> and ip are boson and fermion fields, respectively. 

The fermion mass m/ and trapped frequency w/ are normalized with the boson mass Mb and the 
boson trapped frequency Os, respectively. The coordinates are normalized by £b = (H/Mb^b) 1 ^ 2 ■ The 
coupling constants gs and Hbf are given as 

g B = 4:ira B B £ B \ (5) 
h B F = 2na B F S, B 1 ( 1 + m j 1 )^ ( 6 ) 

where asB and a^i? are the scattering lengths between two bosons and between the boson and the 
fermion, respectively. 

In this work we only consider only zero-temperature, so that we can approximate the system with only 
the condensed bosons and degenerate fermions. Namely the total wave function is written as 

N B 

$(T)={n^)}*/, (7) 

»=1 

where t is the time normalized by fi^ 1 , 4> c is a wave function of the condensed boson, and is a Slater 
determinant of fermions with single particle wave functions ip n . 

The time evolution of the wave functions arc obtained from the variational condition that 



5 ( 2 cLtL{t) =8 ( 2 dr < $( T )|{i-£- - H}\&(t) >= 0. 
J Tl Jn OT 



(8) 



From this condition we derive coupled equations of the TDGP and TDHF equations as follows. 

i^ c {x,T) = |-iv 2 + Er fl (aOj <t> c {x,r), (9) 

i^ n (x,T) = |--i-V 2 + [/ i ,(a ; )| i> n {x,r) (10) 

with 

U B {x) = \<s? +9bPb(x) + h BF p F (x), (11) 

U F (x) = ^m f uj 2 f x 2 + h BF p B (x), (12) 

where r is the time normalized by fi^ 1 and p B and p F are boson and fermion densities which are given 
as 

p B (x) = N B \4> c {x)\\ (13) 

occ 

PF (x) = ]T|^0r)| 2 . (14) 

n 

The number of fermions are too large to solve the above TDHF equations directly, so instead one uses 
the semi-classical approach. In the semi-classical limit (h — > 0) the TDHF equation is equivalent to the 
following Vlasov equation 

±f( x ,p-r) = !^- + ^-V x -[V x U F (x)][V p ])f(x,p;T)=0, (15) 
cIt ( or rrif J 

where f(x,p;r) is the fermion phase-space distribution function defined as 

f(x,p,r)= J d 3 z<^j(x + ^z,r)^(x-^z,T)\^>e- l P z ^\ (16) 

As an actual numerical method we introduce the collective coordinate method and the test particle 
method to solve the TDGP equation © and the Vlasov equation (|15|l . respectively. 
The wave function of the condensed boson <f> c is expanded with the s-wave harmonic oscillator wave 
function u n (x) as 

N baBC -l 

<j> c {x,r) = £ A n e w "u n {bx 2 )e~i» x \ (17) 



where Nb aS e is the number of the harmonic oscillator bases, and b, v, A n and 9 n are the time-dependent 
variables. 

As for the fermion the phase-space distribution function is described as 

f(x,p,r)= [ -^L 5{x- Xl {r)}5{p~p t (r)} 1 (18) 

iV T f=1 



where Nt is number of test-particles per fermion. 

Then we define the Lagrangian with these coordinates as 

L(A n ,e n ,b,u;xi,p i )=N B J d 3 xcj>* c {x)i^-<j> c {x) + J- ^{Pi^i} ~ E T , (19) 

^ i— 1 

with the total energy Et written as 

E T =< >= E B +E f + E B f , (20) 



and 



e b = J a[^;(x){-v + .uw + |w(*.w) ! ] , (2D 

Ebf = -~-h B Fy^PB{xi). (23) 
In the above equations we take A n , Qn^b^v and Xi,p i as time-dependent variables. 

The time evolution of these variables are given by the Euler-Lagrange equations with respect to the 
time-dependent variables. For example, the Euler-Lagrange equations for the fermions with respect to 
Xi and Pj, become the following classical eq. of motion for the test particles, 

4U(r) - (24) 
= -VI^(aj). (25) 

We can also obtain the same equations by substituting ca. (|18J) into ca. (|15|l . Note here that the test- 
particle motion describes the time evolution of the fermion gas, and does not correspond to actual single 
particle motion. 



III. RESULTS 

In this section we show the results of our calculation on the monopole oscillation in the boson-fcrmion 
mixture. In order to know the characteristic behavior of the boson-fermion mixed system, we consider 
only the system where the boson number Nb is much larger than the fermion number Nf (Nb 3> Nf), 
which optimizes the overlap of bosons and fermions. This is because the condensed boson density tends 
to be located in the central region while the fermion density distributes over a large region due to Pauli 
blocking. 



G 

In this work we deal with the system 39 K - 40 K, where the number of the bosons ( 39 K) and the fermions 
( 40 K) are iV& = 100000 and Nf = 1000. The mass difference between the two atoms is omitted, m/ = 1. 
The trapped frequencies are taken as f2g = 100 (Hz) and ujf = 1. The boson-boson interaction parameter 
gs is taken to be gs — 1-34 x 10 -2 , which corresponds to clbb = 4.22(nm) 27]]; in addition we vary the 
boson- fermion interaction parameter Hbf which is measured in units of ho = 7.82 x 10 -3 corresponding 
to cibf — 2.51 (nm) [27j. In the numerical calculation we take the number of the harmonic oscillator 
bases for the condensed boson to be Nb ase — 11 and the number of the test-particle per fermion to be 
N = 100. 



A. Ground State 

Before performing the numerical calculation of collective motion, we explain the method to construct the 
ground state under the above formalism. The ground state is defined as the stationary state with the 
lowest energy. 

In the stationary condition the wave function of the condensed boson satisfies v — 6 n — 0, and dET/db = 
dEx/dAi = 0. We can obtain the wave function of the condensed bosons by repeating the processes 
of solving these stationary conditions and diagonalizing the Hamiltonian matrix with the base of the 
harmonic oscillator wave functions. 

The fermion phase-space distribution function is given by the Thomas-Fermi approximation, 

f(x,p) = 6\nf-e(x,p)] (26) 

with 

e{x,p) = ^—p 2 + U F {x). (27) 
2m/ 

Here we iterate solving the boson wave function and varying the fermi energy fif to give the correct 
fermion number. 

For numerical simulations we distribute the test-particles to reproduce the phase-space-distribution solved 
in the Thomas-Fermi approximation. In this process we have to insure the stability of the ground state. 
In Fig. we show the density distributions of the bosons and the fermions (multiplied one hundred) 
with the long-dashed and dashed lines, respectively. We take the boson- fermion coupling to be Hbf = ^o- 
We see a large overlap region between boson and fermion densities. In addition we show test-particle 
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distribution with Nt — 100 in Fig.^3, where the histogram indicates the test-particle distribution while 
the dashed line denotes the density in the Thomas-Fermi approximation. We see that the two results 
very nicely agree. 

As a next step we examine the stability of the ground state by performing a time evolution starting from 
above the ground state. In Fig. [2 we show the time dependence of the root-mean-square-radius (RMSR) 
for the bosons (Rb) and the fermions (Rf) with long-dashed and dashed lines, respectively. These results 
show that their RMSRs vary only within 0.01 % over the entire time-evolution, and confirm that the 
ground state is very stable. Hence we can get very stable and reasonable ground states for the dynamical 
calculations. 



B. Monopole Oscillations 

In this subsection we show our calculational results for the monopole oscillations. We take the boson- 
fermion coupling to be Kbf — m all calculations in this subsection. 

At the beginning of the numerical simulations we scale the boson wave function <f> c and the fermion phase- 
space distribution function from those in the ground state as follows. By using the scaling parameter Sb 
and s/, we scale the boson parameter b and the position (xi) and momentum coordinates of the i-th 
fermions as 

6(0) = 6 (9) /^i (28) 
Xi(0)=8 f xW, Pi (0)= P ^/ Sf , (29) 

where superscripts (g) represent the coordinates in the ground state. Then we perform the numerical 
calculation of the monopole oscillations with various initial conditions. 
In order to see appearance of the monopole oscillation, we define 

Ax i (r)=R i (r)/R^-l (i = B,F), (30) 

where i?s(r) and Rf{t) are the RMSRs of the bosons and fermions at time r, and R? B and R? F are the 
RMSRs in the ground state obtained by the Thomas-Fermi approximation. 

In Fig. |3 we show the time-dependence of Ax B (a) and &.xp (b) with the initial conditions Axb{Q) = 0.1 
(sb = 1.1) and Axf(0) = 0.1 (sf = 1.1) (in-phase). In addition we exhibit results with other initial 
conditions: Ax B {0) = 0.1 (s b = 1.1) and Ax F {0) = -0.1 {s f = 0.9) (out-of-phase) in Fig.^b, Ax B (0) = 
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0.1 (sf, = 1.1) and Ax^(0) = (s/ = 1.0) in Fig. Eli- The condensed bosons oscillate independently of 
the fermions, so that the difference between Axb at these three initial conditions is not visible. 
In the above three cases we see that the fermion oscillations are forced vibrations with a beat, which 
gradually becomes blurred. This blurring cannot be explained within the RPA. 

Before discussing this blurring we calculate with another initial condition: Axb(0) = (s& = 1.0) and 
Axp(0) = 0.1 (sf = 1.1) . With this initial condition the boson oscillation is weak and does not strongly 
affect the fermion oscillation, see Fig. 0] The boson oscillation has a beat which appears in early time 
and blurs in latter time stages, and the fermion oscillation has a strong damping. In Fig. |2t we show 
the results of the calculation at the same initial condition, but with the boson motion frozen; namely the 
fermions move in the fixed potential Up(x) in the ground state. The damping behavior of the fermion 
oscillation is almost the same as that in Fig. 0Jd, except that the fermion oscillation is monotonously 
damped while the oscillation in Fig. 0Jd undergoes one beat. 

This damping process causes the blurring beat in the fermion oscillation shown in Fig. [3J Hence this 
damping process must play a significant role in collective motions in the boson-fermion mixtures. Next we 
inspect the origin of the damping by focusing on monopole motions with the initial condition Axb(Q) = 0. 
Here we examine the time-evolution of the energy levels. In Fig. [3] we show the time-dependence of Eb 
(long-dashed), Ef+Ebf (dashed) and Ebf (solid), which are defined in eqs.J2UH2III- Note that Ef+Ebf 
corresponds to the sum of the fermion single particle energy. 

First we note that the boson-fermion interaction part, Ebf, is much smaller than the boson part, Eb, and 
the fermion part, Ep- Second the boson part Eb and the fermion part, Ep + Ebf, do not significantly 
vary in the time-evolution process. After the damping (r > 100) the energies, normalized by MIb, are 
given as 

E B = 7.04 x 10 5 (7.04 x 10 5 ), (31) 
Ep + Ebf = 1-57 x 10 4 (1.54 x 10 4 ), (32) 
E BF = 1-54 - 1.61 x 10 3 (1.57 x 10 3 ), (33) 

where the values in the brackets indicate the energy contribution to the ground state. The excited energy 
of the fermions is about 300 (MIb) while the energy transfer is about 70 (MIb) through the boson-fermion 
interaction. This fact implies that the damping is not caused by the fermion energy loss, and in addition 
that this oscillation state is a many-particle many- hole state because the excited energy of the one-particle 



and one- hole states is about 2M7# in the monopole oscillations. 

Further we calculate the spectrum which is obtained by the Fourier transformation 



F(co) 



[ ' dTR 2 (r)e WT 



(34) 



In order discuss the oscillations, we vary Ti and T/. 

In Fig. we show the results of the spectra obtained by integrating over the regions < r < 80 (damping 
period) for the bosons (a) and fermions (b), and over 100 < r < 200 (after damping ) for the bosons (c) 
and fermions (d). In addition we plot the spectra for fermion oscillation with the boson motion frozen, 
with the dashed lines in Fig. and Eli. 

In Fig. |HJd we see that the peak position for the fermion spectrum isw~ 1.91, and does not vary in any 
of the time regions while the width becomes slightly narrower after the damping is finished (Fig. Eli) . 
On the other hand there are two peaks for the boson spectra; at the frequency u « 1.91 and 2.23 in 
the fermion damping process. After the damping, there is only one boson oscillation mode at w w 2.23, 
similar to the fermion oscillation. 

We can consider the intrinsic frequencies of the boson and fermion oscillations to be u) b M = 2.23 and 
ui^ M = 1.91 in our calculation. The boson intrinsic frequency agrees with u> = y/5 which has been shown 
in the system only including the condensed bosons |2j. The boson intrinsic mode is almost independent 
of fermions because the boson number is much larger than the fermion number. In addition the frequency 
of the fermion oscillation is the same with or without the boson motion. 

In the short time period both the boson and fermion oscillations can be approximately described within 
the linear response theory as 

Ax B ~ A bb cos(u; b M T) - A bf cos(u f M T), (35) 
Ax F » A fb cos(uj b M T) - A f f cos(oj f M T). (36) 

If the amplitude of the boson oscillation in the initial condition is not small, the boson oscillation negligibly 
involves a mode with the fermion intrinsic frequency, \A bb \ 3> and the fermion oscillation becomes 

the usual forced oscillation \ Af b \ i=s \ Aff\ as shown in Fig. EJ 

When Axb(0) — such as in Fig.^J the fermion oscillation does not have a mode with the boson frequency 
(| Af b \ sa 0). In the beginnings these fermion oscillations trigger the boson oscillations, and then the boson 
oscillations dominantly hold the mode with the fermion frequency, A bb 3> A b f. As the time evolves, the 
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amplitude of the boson oscillation becomes larger, but the amplitude of the fermion oscillation contrarily 
decreases, and the A^f is gradually reduced by the damping of the fermion oscillation. The amplitude of 
the fermion oscillation becomes too small to influence the boson oscillation, the boson oscillation loses 
the mode with the fermion frequency, and its beat finally disappears. 

The behavior of the fermion oscillation in the full calculation is the same as that when the boson motion 
is frozen, and the fermions do not lose energy through oscillation. Hence the fermions move almost 
independently each of other, and the damping of the fermion oscillations must come from properties of 
the fermion single particle motion. 

The boson density is distributed in a smaller region than the fermion density, so the fermion potential 
does not have simple harmonic oscillator shape, which is shown with the solid line in Fig.0 The Thomas- 
Fermi approximation, which is available in the case of the large boson number, gives the boson density 
as ps ~ — x 2 )/2gB, where fiB is the boson chemical potential. Then the fermion potential Uf is 

approximately written as 



Thus the fermion potential is separated into two region, the inside and outside of the boson populated 
region. 

In the single harmonic oscillator potential the fermion motion is always harmonic in both the semi-classical 
and quantum approaches, but the above fermion potential, which is separated into the two regions, leads 
to anharmonicity when fermions move in the both regions. 

In the quantum calculation this damping is explained with the phase factor of the fermion single particle 
wave functions given by their single particle energies. If the boson motion is frozen, the wave- function of 
the n-th fermion is written as 



where Vi is the wave-function of the i-th orbit of the ground state, Si is its single particle energy, and a\ 
is a coefficient. The expectation value of x 2 is 




(37) 




(38) 




(39) 
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with 

Cij = ^ a< n) o< n >* / ^(^(x). (40) 
This cocfBcicnt is usually negligible when \i — j\ > 2, and then 

< a; 2 >w ^ {(7^ + C i , i+1 e i ( £i + 1 - £ *) T . + Q+ M e- , < e «+ 1 - e '> r } . (41) 

i 

When the fermion potential is only the harmonic oscillator potential, e^+i — e,; = 2wy, and the oscillation 
of < x 2 > does not exhibit a damping. As shown before, however, the fermion potential includes two 
kinds of harmonic oscillator potentials with the different frequency. In addition the potential in the 
boson region is not exactly that of the harmonic oscillator. Then the difference of the single particle 
energies, £j+i — has dependence on i, and the fermion monopole oscillation includes various vibration 
modes with different freqencies, and its oscillation shows a damping. It was shown in Ref. |28j that 
an oscillation of population difference in two-component fermi gase has a damping because of a similar 
multimode dephasing. 

In the semi-classical calculation, the frequencies of modes including the oscillation is continuous, though 
they are discrete in the full quantum calculation. In the system with large particle number, however, 
the semi-classical calculation well describes the time evolution of the phase space distribution. In the 
harmonic oscillator potential, especially, this time evolution completely agrees with that in the quantum 
calculation because higher order terms with respect to h 25], which are omitted in ea. l|15|) . exactly vanish. 
Hence our calculations sufficiently describe the actual time-evolution in the mean-field level. 
In numerical simulations such oscillation properties are described by motions of the test-particles. A test 
particle passes over both regions, inside and outside of the boson region. Here we consider oscillation 
behavior of radial distance of this i-th test particles as \xi\. In a non-harmonic oscillator potential the 
period between successive maximum \xi\ depends on the orbit. This anharmonicity gradually reduces 
the amplitude of Ax. 

As mentioned before, the amplitude is about Ax ~ 10~ 4 in the RPA calculation and the excited energy 
is about 2ft£lB- In the present calculation, however, the excited energy when Axf(0) = 0.1 (sf = 1.1) is 
about SOOHUb- Oscillation states such as those in the present calculation are multi-particle and multi- 
hole states. Under an anharmonic potential, all test-particles do not oscillate with the same period, and 
hence this anharmonicity is the cause of the observed damping of Aa; p ■ 
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C. Amplitude and Boson- Fermion Coupling Dependences 

In this subsection we investigate the oscillation behavior by changing the initial amplitude of fermion 
oscillation, Axf(0), and the boson- fermion couplings, hsF- 

First we examine the amplitude dependence of the damping. In Fig. |H1 we show the results with the 
initial conditions that Axb(0) = and Axf(0) = 0.01 (small amplitude). The excited energy of this 
oscillation is about 30 (MIb), which is still much larger than the one-particle one- hole excited energy. 
The strong damping also appears there even with such small amplitude oscillation, which becomes about 
10 -3 after the damping. The test particle method cannot however give sufficient accuracy to exactly 
describe oscillations with such small amplitudes after the damping. Therefore, the more detailed study 
lies beyond the scope of this work. 

In Fig. we show the results with the initial conditions that Axb{0) = and Axp(0) = 1.0 (large 
amplitude); we do not see such strong damping though the excitation energy is about 1800 (hflg). In 
Fig. ^| we show the spectra of the above boson and fermion oscillations for the large amplitude. The 
boson and fermion spectra have peaks at quite different frequencies. The boson spectrum has a large 
peak at the boson intrinsic frequency u rs uj^ m (a), while the fermion one has a peak at uj w 2 (b), 
which is the frequency of the monopole oscillation with only one harmonic oscillator potential. In the 
large amplitude oscillation most of the fermions have large kinetic energy in the central region and pass 
through this region quickly, so that the fermions mostly stay outside the boson density region. In large 
amplitude oscillations, therefore, the oscillation frequency agree with that only in the trapped potential 
and, the rapid damping does not appear. The final amplitude after the damping must be much larger 
than that when Axf(0) =0.1. 

Next we vary the strength of the boson- fermion coupling Hbf, and perform calculations of the oscillation 
at the same initial condition Axb(0) = and Axf(0) = 0.1 . We show the results with hsF/ho = 2 in 
Fig.[TT]and h BF /h Q = -1 in Fig.lTl 

In Fig. ^2 the features of the oscillation which have been seen in Fig. 0] are more clear; the damping of 
the fermion oscillation and the blur of the beat in the boson oscillation occurs sooner with hsF = 2/io 
than hsF = ha. 

If the boson- fermion coupling is negative, the oscillation behavior is a little different. When hsF = — ho 
in Fig. 1121 the fermion oscillation undergoes beating even after the motion is damped. Additionally 
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the amplitude of the boson oscillation is much larger than that of hsF > 0. This beat of the fermion 
oscillation after the damping is caused by the rather large amplitude of the boson oscillation. 
As the repulsive interaction between bosons and fermions is increased, the damping of the fermion oscil- 
lation is faster, but the amplitude of the boson oscillations do not increase. In the attractive interaction, 
however, the damping is still rapid, but the amplitude of the boson oscillation is large. We believe this 
difference comes from the overlap region between boson and fermion densities. 

When the fermions feel a repulsive force in the boson density region (see Fig. [7J) , the fermions stay outside 
of this boson region and have little influence on the boson oscillation. Then the amplitude of the boson 
oscillation does not become larger, and its beat becomes blurred soon in spite of the stronger interaction. 
When the interaction is attractive, on the other hand, many fermions moves inside of the boson density 
region, and then boson and fermion motions are more significantly coupled. 

IV. SUMMARY 

In this paper we study the collective monopole oscillation of bose-fermi mixtures, where the number of 
bosons is much larger than that of fermions, by solving directly the time-dependent equations. When 
the initial amplitude is about 10 % of the ground state RMSR, we find a rapid damping of the fermion 
oscillation and a beat in the bosonic vibration at the early time stage while this damping and beat almost 
disappear in the later time stage. This rapid damping does not exist at zero temperature in the system 
including only one kind of boson 29] ; it is a typical feature in the bose-fermi mixtures. 
Our analysis shows this damping is caused by the anharmonicity of the multi-particle and multi-hole 
states. Since the fermion mean-field is separated into two parts, inside and outside of the boson density 
region, fermions feel very different forces in these two regions and their motions are not harmonic. Then 
these motions becomes chaotic, and the oscillation amplitude decreases. 

In this work we focus on the oscillation with the initial condition Axb(0) = 0, where only the fermions 
have initial motion at the beginning. In this case we can clearly see damping properties of the fermion 
oscillation because the boson oscillation has very small amplitude. Even if the amplitude of the initial 
boson oscillation is larger, the damping also plays a significant role. When the fermion oscillations have 
a large amplitude, it always tends to cause damping, but the amplitude after the damping is enhanced, 
again, by the boson oscillation. The fermion oscillation repeats the damping and enhancement, and its 
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beat gradually blurs; the oscillation behavior is not recursive. 

This damping process cannot be described with RPA, where the motion is assumed to be harmonic. As 
shown in Ref. 15| the excitation energy obtained in RPA corresponds to oscillations with the amplitude 
of 10~ 4 of RMSR. In actual experiments, however, the amplitude is not so small, the oscillation states 
are multi-particle and multi-hole states, and hence we must solve the time-dependent process directly. 
Even if the boson-fermion interactions are attractive, both the boson and fermion oscillations have an 
oscillation behavior similar to those in the repulsive interaction in the damping time stage. In the 
attractive interaction, however, the amplitude of the boson oscillation slightly increases and the fermion 
oscillation has a small beat after the damping time stage. 

In this work we do not show calculations with stronger attractive force between bosons and fermions. If 
the attractive force is stronger, the frequency of the fermion oscillation becomes larger, and finally agrees 
with that of the boson intrinsic oscillation; where the two oscillation must make a resonance. We have 
investigated this phenomena, and found resonances which have quite different behaviors 
Furthermore we do not take into account two body collisions and thermal bosons. Since two fermions 
do not directly interact, dissipation of fermions occurs through collisions between fermions and thermal 
bosons. In the system Nf, 3> Nf at zero temperature, however, the thermal bosons are very few, and then 
the damping process shown in this work must be dominant in the early time stage. After this damping 
the collisional damping may play a role in the time-evolution process. In future we would like to introduce 
this collisional process in our simulation. 
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FIG. 1: The density-distribution of ground state of the boson (K 39 ) and fermion (K 40 ) mixed system with 
Iibf = ho- The results of Nb = 10000 and Nf = 1000 arc shown in Fig. la. The dimensionless distacnce x is in 
units of £b. The long-dashed and dashed lines indicate the densities of boson and fermion, respectively. In lower 
column the histogram show the results of the test particle. 
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FIG. 2: The time-dependence of the root-mean-square radii normarized by £b for Nb = 100000 and Nf = 1000 
with Hbf = ho. The dimensionless time r is in units of O^ 1 . The long-dashed and dashed lines indicate the 



boson and fermion root-mean-square radii (RMSR), respectively. 
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FIG. 3: Time evolution of Axb (a) and Axf (b) with Iibf = ho at the initial condition that Aa^s(O) = 0.1 
and Axf(0) — 0.1. Time evolution of Axf at the initial conditions that Sb = 1.1 and Axf(0) = —0.1 (c), and 
Ax F {0) = (d). 
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FIG. 4: Time evolution of Ax for boson (a) and fermion (b) with Hbf = ho at the initial condition that Axb(0) = 
and A^f(0) =0.1. In the last column time evolution of Ax for fermion (c) with freezing the boson motion at the 
same initial condition. 
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FIG. 7: Fermion potential at the ground state. The dashed, solid, and chain-dotted lines show those at Kbf = 2/io, 
Iibf = ho and Hbf = — ho, respectively. The dotted line indicates the harmonic oscillator potential. 
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FIG. 8: Time evolution of Ax for boson (a) and fermion (b) with \ibf = ho at the initial condition that Aa^s(0) = 
and Aa: F (0) = 0.01. 




FIG. 9: Time evolution of Ax for boson (a) and fermion (b) with Hbf = ho at the initial condition that A^s(O) = 
and Aa: F (0) = 1.0. 
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FIG. 10: Spectra of the boson oscillation in upper columns (a) and fermion oscillation (b) at the initial condition 
that Aa;s(0) = and Axf(0) = 1.0. The spectra are deduced from the evolution in < r < 200 in the oscillation 
shown in Fig.[§] 
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